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■ Abstract. We implement complex scaling of Faddeev equations using hyper- 

spheric coordinates and adiabatic expansion. Complex scaling of coordinates 



o 



I allows convenient calculations of three-body resonances. We derive the nec- 

2 ' essary equations and investigate the adiabatic spectrum at large distances. 

We illustrate the viability of the implementation by calculations of several 
three-body resonances: a 0"*" resonance in a model benchmark system of three 
identical bosons; the 2"*" resonance in the ^He nucleus within the a + n + n 
^ ■ model; and the two O"*" resonances in ^■^C within the three-a model. 

. ^ 1 Introduction 

^ , Although the quantum-mechanical three-body problem has a long history 1 it 

is still generally considered as "difficult" . Recently an effective symbiotic method 
was introduced [21 where the Faddeev equations |31 are solved using the hyper- 
spheric adiabatic expansion . The method allows analytic asymptotic solutions 
and thus is especially suitable to describe weakly bound and spatially extended 
systems. During the last years it has been successfully employed to investigate 
ground state properties and fragmentation reactions of three-body halo nuclei 
^He and ^^Li ^MU\ 

Investigation of three-body systems in basically all fields of physics led very 
soon to the question of the three-body continuum spectrum. Investigations of 
three-body resonances is today a hot topic in both experimental and theoretical 
physics. 

In the early 70's the complex scaling method was introduced |H1 IH] to cal- 
culate resonances, where the radial coordinate, and therefore the Hamiltonian, 
is rotated into the complex plane. This transformation does not change the en- 
ergies of the bound states, while resonances appear as additional bound states 
with complex energy values. Thus the usual numerical techniques used to com- 
pute bound states can also be employed to calculate resonances ^1 ^2 El Ej 
(see also El El E] for recent applications to three-body systems) . 
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A popular alternative coordinate space method to calculate resonances is the 
complex energy method where the poles of the S'-matrix are searched for. This 
method requires matching of the calculated wave function with a certain large 
distance asymptotics, typically a Bessel function, but sometimes a more com- 
plicated function. Furthermore the asymptotic functions diverge exponentially, 
and an accurate calculation at large distances becomes soon a delicate numerical 
task, especially for the long-range Coulomb potentials [T8] . 

The great advantage of the complex scaling method is that it is possible 
to avoid matching at large distances and to employ instead a simpler boundary 
condition of vanishing wave function or, alternatively, a basis of square integrable 
functions. 

In this paper we implement the complex scaling method within the hyper- 
spheric adiabatic approach to the Faddeev equations. We first derive the neces- 
sary equations and investigate the asymptotic adiabatic spectrum. We then il- 
lustrate the viability of the method by applications to three systems: i) a model 
benchmark system of three bosons; ii) the halo nucleus ^He within the three- 
body a + n + n model; and iii) the three-o system where the long range Coulomb 
interaction is of importance. 



2 Complex scaling of hyper-spheric coordinates 

In this chapter we reintroduce some pertinent definitions from the complex scal- 
ing and hyper-spheric adiabatic methods, perform the scaling of hyper-spheric 
coordinates, derive the formulae necessary for the numerical implementation of 
the method, and discuss the large-distance asymptotics of the hyper-angular 
adiabatic spectrum. 



2.1 Complex scaling in a two-hody system 

A resonance or, in the stationary scattering language, a decaying meta-stable 
state corresponds to a pole kQ = q — oi the S-matrix located in the IV-th 
quadrant of the complex momentum. The corresponding radial wave-function 
fko asymptotically grows exponentially 

(X e^^^e^r (1) 

The complex scaling method is designed to avoid this exponential growth by 
scaling the coordinate r with a complex exponent e*^ 

r ^ re'^. (2) 

This changes the asymptotics of the radial wave-function into 

fk OC Q^\^Vi\r cos(e - ip) ^-\ko\rsin{e~ if) 

where —ip is the argument of the complex momentum (/cq = l^ole"**^)- If 
d > i-p this wave-function is exponentially vanishing and the resonance effectively 
becomes a bound state and thus numerical methods designed for bound state 
problems could also be employed for resonances. 
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The exponentially vanishing wave-function Q allows a simple bounding box 
boundary condition 

/(w) = 0, (4) 

provided r^ax is sufficiently large. The continuum spectrum solution for the 
complex scaled radial function at large distances is 

f{r) ocsm{kre'^ + 6), (5) 

where 5 is the phase shift. The boundary condition with this function then 
leads to a trade-mark rotated spectrum of the continuum states lying along a 
line rotated by —9 {—20 in the energy plane) 



.„_^e-^ n = l,2,... (6) 

^max 

Although the wave-function © is exponentially vanishing, the fall-off rate is 
controlled by the parameter sin(0 — (/?) and therefore r^^^ must satisfy 

|A;o|rmax sin(6' - 93) > 1, (7) 

and thus is typically larger than for bound states, where only 

|A:o|rmax » 1 (8) 

is needed. Also for exponential, exp(— r/6), and gaussian, exp(— r^/6^), potentials 
the range h of the potential is increased after complex scaling as, correspondingly, 
6/ COS0 and cos 26. For gaussian potentials the scaling angle therefore can 
not exceed 9 = 7r/4. 

This increase of the range of the scaled potential is counterproductive as it 
necessitates numerical calculations to larger distances. 

2.2 Hyper-spheric adiabatic approach 

We first define the usual Jacobi coordinates Xj and in terms of masses rui and 
coordinates rj of the three particles (i=l,2,3) 

Xj = ^/JH{rj - Ffc) , (9) 

/ m,r,' -|- m^rfc \ 

= V^i^-- m,+m, ) ' 

1 mj-nik 1 mi{mj + rrik) 

m nij + rrik mrrii + mj + rn^ 

where {^, k} is a cyclic permutation of {1,2,3} and m is an arbitrary mass. The 
hyper-radius p and the hyper-angles Ui are now defined as 

psin(ai) = Xi , pcos(Qi) = m . (11) 

A set of hyper-spheric coordinates consists of the hyper-radius p and a set of five 
angles f2i - the hyper-angle Oj and the four directional angles of Xj and y^. Here 
i distinguishes different sets of Jacobi coordinates. 
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The kinetic energy operator T is then given as 

T = Tp + 



2m 



1 15 



A' 



sin(2ai) daf 



sin(2Q:j) - 4 + 



+ 



ll 



sin^(aj) cos2(aj) 



(12) 
(13) 
(14) 



where la;, and \y^ are the angular momentum operators related to Xj and yj. 

The hyper-spheric adiabatic method allows effective solution of the three- 
body problem with the Hamiltonian 



2mp^ 



(15) 



where Vi is the two-body potential between particles j and k. The coordinates are 
divided into the "slow" hyper-radius p and "fast" hyper-angles Q (here we refer 
to any of the three sets of angles). For every fixed hyper-radius the eigenvalue 
problem 

.2 3 



(16) 



is solved and the eigenvalues Xn{p) and eigen- functions ^niPi ^) £^re calculated 
as function of p. The total wave-function is then expanded in terms of these 
eigen-functions 



p5/2 

The expansion coefficients fn{p) satisfy a coupled system of differential equations 
where the angular eigenvalues serve as effective potentials 



^2 A„ + 15/4 2mE 



dp' 



+ 



p^ 



d 
dp 



fn', (18) 



where E is the three-body energy and P and Q are the non-adiabatic terms 
defined as 



Qnn'ip) = J ^liip, ^)-g^^nf{p, i^)df2, 

where f denotes hermitian conjugation. 
2.3 Scaling of hyper-spheric coordinates 

We now perform the complex scaling of the Jacobi coordinates 



(19) 



(20) 
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The hyper-radius p is then scaled with the same angle 9 while the hyper-angles 
Q remain real as they only depend upon ratios of the Jacobi coordinates 



pe 



(21) 



We shall in the following always maintain real p writing the complex rotated 
variable explicitly as pe*^. The angular eigenvalue problem ()16() acquires then 
complex potentials 



2ie 



2mp 



i=l 



(22) 



and therefore the eigenvalues A„ and eigen-functions ^„ also become complex. 

We solve the complex-scaled eigenvalue problem H22() by using the hyper- 
spheric harmonics as a basis set 0. Since the hyper-angles are not affected by 
the complex scaling these functions remain real. The operator TC at the left hand 
side of the eigenvalue equation , 



i=l 



(23) 



is in this basis a symmetric complex matrix, Ti'^ = TL (here denotes transpo- 
sition). A non-scaled hamiltonian, ^ = 0, is a real symmetric matrix, that is a 
hermitian matrix. 

A complex matrix Ti has in general different left, &iid right, eigenvectors 
corresponding to a given eigenvalue A 



(24) 



The left Xn and right eigenvectors corresponding to different eigenvalues A,^ 
and A„/ are orthogonal in the following way 



(25) 



For a symmetric matrix there is an apparent connection between the left and 
right eigenvalues 



X 



(26) 



Therefore for the complex scaled angular eigen-vectors the scalar product rule 
takes a somewhat unusual form 



(27) 



and thus the hermitian conjugation in the definitions of P and Q in must 
be substituted by transposition in the case of complex scaling. 

The hyper-radial equations (fT8|) remain largely unchanged, but, of course, the 
eigenvalues A„ together with the non-adiabatic terms P and Q are now complex. 



6 Complex scaling of the hyper-spheric coordinates... 

Similar to the two-body three-body resonance corresponds to a pole 

of the three-body 5'-matrix at a momentum k = q — i-j {k = y^2rnE/n?). At this 
point all hyper-radial functions fn have the diverging asymptotics 

fn oc e^'^Pe'TP. (28) 

The complex scaling of the hyper-radius leads in complete analogy to the 
two-body case to an oscillating and exponentially vanishing resonance function 

and again the boundary condition 

fniPmax = 0) (30) 

leads to the rotated discretized continuum states 

TTj 



j = l,2,.... (31) 



2.4 Faddeev equations. 



The angular wave-function ^ is now represented as a sum of three components 
(p = -\- (^(2) _j_ ^(3) -^i^ic}! satisfy the three Faddeev equations 

where k=l,2,3. One can easily verify that the sum of these three Faddeev equa- 
tions is equivalent to the Schrodinger equation (|22j) . 

The Faddeev equations H32|) can be conveniently written in a matrix form as 

Fcj) = X(j), (33) 

where (j) is a column vector of the three Faddeev components 

(34) 

V 0^') J 

and the Faddeev operator F is equal 

F = T + VR, (35) 
where T is the kinetic energy matrix 

/ A"^ 

T = I , (36) 

\ yl2 

V is the potential matrix 

V = I V2 I . (37) 




Vi 






















V3 
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and R is the projection matrix 

i? = I 1 1 1 I . (38) 




If we again use hyper-spheric harmonics as a basis set then each of the ele- 
ments in the above matrices becomes in turn a matrix in this basis. 

The matrix R acting on the Faddeev vector H34|) apparently produces a vec- 
tor with the total function <P as the three components. The matrix R has the 
following properties 

R"^ = 3R, (39) 
R^ = R 

The scalar product rule (|27(1 for the Faddeev functions becomes 

cgR^RK' = 5nn' ^ (40) 

T 1 
(j)nR4'n' = gl^nn'- 

The factor 1/3 can be hidden into the normalization by redefining the wave- 
functions such that 

4>lR(j)n' = Snn'- (41) 

This product rule should be used in the definitions H19() of the non-adiabatic 
terms P and Q if one uses the Faddeev three-component functions: 

Pnn'ip) = I cgR^dQ, (42) 



Qnn'{p) = l&^dn. 



Actually, the product rule with the i?-matrix (|4H) represents the ordinary 
product rule for the Faddeev operator F. Although the Faddeev matrix is not 
hermitian, it has the important symmetry property 

RF = F^R, (43) 

which leads to a simple connection between the left, and right, (f), eigenvectors 
of the matrix F 

(/?^ = (t>'^R, (44) 

which immediately leads to H41|) . 

Equation (|44j) can be proved by multiplying i|33j) with R from the left, trans- 
posing it and using the relation H43() . Apparently, since we did not assume that 
the matrix V is real, the complex scaling does not affect the validity of the above 
expressions. 

Note that for the non-scaled real F the property (|43)) means that the eigen- 
values of F are real although the matrix itself is not hermitian. 
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2.5 Phase of the angular eigen-functions 



The complex eigen-functions ^„ obtained from the complex scaled eigenvalue 
equation (|^^ are as usual defined up to a phase factor. We can therefore multiply 
the eigen- function <Pn by an arbitrary phase function e*-^"^''\ where Xn is real. 
This changes the non-adiabatic terms P and Q in 1)19^ as 



Pn 
Qn 



Pnn' ~l~ '^Xn'^nn' 7 

Qnn' + '^ix'n'Pnn' + {ix'n' " Xn') ^n 



(45) 
(46) 



(where primes denote differentiation with respect to p) and correspondingly, via 
equation ((TH|l . leads to different radial functions /„. However, the total wave- 
function 'F (X ^nfn in l(T7|) does not change as the new set of radial functions 
merely corresponds to a trivial phase shift 



fn 



fn 



(47) 



such that the product of the angular and the radial function remains unchanged 



^nfn ^ e'^'^^ne-'''"^''^fn = ^nfn- 

Indeed the right hand side of the radial equations (fTH|) transforms as 



(48) 



d 

'^Pnn' I" Qr 

op 



-«X„' ( 



2P„ 



d_ 
dp 



~l~ Qn r 



fn 

d 



fn' + ix'n' + ^Xn') ^nn' fn' + 2iXn'Snn' ^fn') 



(49) 



The last two terms, however, precisely cancel with the corresponding terms from 
the second derivative at the left hand side of the radial equations (|18p 



' fn 



9p2 



o2 o 
fn = - ^X'nfn " Xn fn " ^iXng^fn)- (50) 



Thus the phase of the hyper-angular eigen-functions can be chosen arbitrarily 
as the total wave-function is independent of this phase. 

In the non-scaled case the diagonal elements Pnn are identically equal zero. It 
is convenient to chose the phase such that this property holds also for the scaled 
case. One can easily verify from H45|) that the phase 



Xnip) = exp ( i / Pnn[p')dp ] , 







(51) 



precisely leads to a vanishing diagonal P as in the non-complex-scaled case. 

The freedom in choosing the phase can also be exploited in order to eliminate 
a given non-diagonal term Qnn' ■ The phase should in this case instead be chosen 

as 

' i 
2 



exp 







P-^,{p')Qnn'{p')dp' 



(52) 



as seen from 
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2.6 Asymptotics of the eigenvalues 

In the angular eigenvalue problem ((22j the two-body potentials are multiplied 
by p^. At the origin, p=0, they therefore vanish from the equations and the 
eigenvalues then start from the free spectrum K{K + 4), where K is even or 
odd integer. The next term in the p expansion is proportional to p"^ as in the 
non-scaled case j2j. 

The large distance asymptotics can be investigated by use of zero-range po- 
tentials. For our purpose it is also sufficient to assume that only one pair of par- 
ticles interact with each other. This simple model will show all different types 
of asymptotics except for the Efimov effect |19| . However the latter is not the 
subject of the current investigation. 

In this case the physical (non-spurious (2j ) solution to the Faddeev equations 
(|32jl contains only one non-zero Faddeev component in (|34|1 . namely the one 
where the two-body potential acts. 

The zero-range potential is vanishing everywhere except for the origin. The 
physical (one-component) solution to the equation H32|l is then a simple function 

vr 

(j) = sin(i/(a - -)), (53) 

where u = \/X -\- A. The zero range potential enters as a boundary condition at 
the origin [20] 

_^cos(^ (54) 

where p is the reduced mass of two particles in units of m, and a is the two-body 
scattering length ^ . 

For p ^ \a\ this equation has several types of solutions. First there is a series 
of asymptotic solutions where sin(z/^) is zero 

z.„ = 2n, A„ = (2n)2-4, n = l,2,... (55) 

which is the free hyper-spheric spectrum K{K -|- 4) with K = In — 2. 

Another type of solution may exist when cot(i/^) remains finite but the 
complex quantity v instead is proportional to p. Suppose that the scattering 
length is negative, a < 0, corresponding to a bound two-body state with binding 
energy B = h^/{2mpLa^). In this case there exists an asymptotic solution to 

= J_^e^(e+-/2), (56) 

since in this case cot(i^^) — > i. 

Therefore a bound two-body state always results in an asymptotically 
quadratically diverging eigenvalue X = v"^ — A 

^^_2^^2^2.,^ (57) 

For positive scattering length, a > and 9 < 7^/2 only the free asymptotic 
solutions to the equation equation (jSl]) exist. Positive scattering length 

^The sign convention for the scattering length is kcoi(5) = 1/a + O(fc^) 
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corresponds to a resonance (more precisely a virtual state) which is located at the 
complex momentum k = —i/a, that is with the argument — 7r/2 in the complex 
momentum plane. The energy of the virtual state is given 

When a > and the scaling angle 6 is larger than 7r/2 yet additional asymp- 
totic solution to (|5l|) appears 

^ = ^^S'~^m, (58) 

That is when the scaling angle 9 is larger than the angle of a resonance (7r/2 in 
the case of a virtual state) the latter effectively becomes a bound state resulting 
in an angular eigenvalue similar to ()57() but with the resonance (or virtual state) 
energy Eq instead of the bound state energy 

A = ^p^e^^^ (59) 

Thus the asymptotic spectrum of hyper-angular eigenvalues primarily con- 
sists of the free spectrum In addition each two-body bound state gives a 
quadratically diverging eigenvalue ((F7|) , and finally each two-body resonance (or 
virtual state) whose argument is less the scaling angle 9 also provides a quadrat- 
ically divergent eigenvalue H59() . 

The eigenvalues enter the hyper-radial equations (|18|) as effective potentials 
A/yO^. After divided by the free eigenvalues converge to the three-body 
threshold at zero energy. These eigenvalues asymptotically describe the gen- 
uine three-body states where all three particles are away from each other. The 
eigenvalues (fFTj) corresponding to the two-body bound state after division by 
converge to a lower threshold at the energy of the two-body bound state. These 
eigenvalues asymptotically describe a configuration where two of the particles 
are in the two-body bound state while the third particle is in the continuum. 
Finally, the two-body resonances and virtual states which after complex scaling 
appear as bound states produce the same type of asymptotic solutions as the 
true bound states. 



3 Numerical examples 

In this chapter we investigate several resonances in different three-body systems 
using the formulated complex-scaled hyper-spheric method and compare with 
other methods. The first application is a benchmark calculation of a resonance 
in a model three-boson system |15| with short-range potentials. The second is 
a narrow 2+ resonance in the halo nucleus ^He. The third application is the 
two lowest 0^ resonances in a system of three a-particles where the long-range 
Coulomb forces are particularly important. 

First we solve the angular eigenvalue problem (|22)) by expansion in hyper- 
spheric harmonics and obtain thereby several of the lowest eigenvalues A„(p) 
as functions of p. For the system of radial equations (|18|) we employ the finite 
difference method. A mesh of about 100 points in hyperradius from zero to 
a certain maximal distance Pmax 

is chosen such that the density of points is 
highest in the regions were the non-adiabatic terms P and Q contribute most. 
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Figure 1. The three lowest complex angular eigenvalues for the model system of three identical 
bosons with scaling angle ^=0.35. 

We then approximate the second and first derivatives in H18|) by a five point finite 
difference approximation assuming that the radial functions are zero at the origin 
and at pmax ■ The system of differential equations (|18() then turns into a complex- 
matrix diagonalization problem. Numerical diagonalization of the matrix gives 
the complex eigenvalues which include all bound states, resonances and also 
discretized continuum states. The resulting eigenvectors give the values of radial 
functions at the mesh points. 

The shooting method with direct Runge-Kutta integration of the radial equa- 
tions gives a better precision but requires a separate search for each of the com- 
plex eigenvalues. Thus the preferred strategy could be first to get an overal 
picture of the complex spectrum by the finite difference method and then to 
calculate selected eigenvlaues with higher precision by the shooting method. 

For each of the following examples we show the angular eigenvalues, the 
complex energy spectrum and the hyper-radial functions of the resonances. As 
these calculation only serve the purpose of illustration we have not optimized 
the hyper-radial mesh to achieve the maximum accuracy. With about 100 hyper- 
radial points the accuracy is about 10 keV in the complex energy. 

3.1 A model three-body system 

In reference a resonance has been calculated in a model three-body system 
of identical bosons with the nucleon mass. The model potential acts only on the 
s-waves and consists of a gaussian attractive pocket and a barrier 



V{r) = -55MeVexp(-0.2fm"V2) + l.SMeV exp(-0.01fm-2(r - 5fmf). 



A three-body resonance was found in |15j ai E = —5.9525 — 0.4034z MeV by 
solving two-dimensional equations. 
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Figure 2. The complex energy spectrum of the model three-boson system with scaling angle 
6^0.4. The bound state at E = -37.22 MeV is not shown. The resonance at _E = -5.96 - 
0.40i MeV is indicated by an arrow. The two lines, Im(_E) = — tan(26')Re(i7) and lm{E) = 
-tan(26')(Re(S) - E'^^^) (where -E^^' = -6.76 MeV is the two-body bound state energy), 
represent the exact rotated continuum spectrum. 



In our calculation we adopt the nucleon mass m=939 MeV as the mass unit 
and use hc=197.3 MeV fm. The three lowest angular eigenvalues obtained by 
solving equation (|22j) for the complex scaling angle 9=0.4: are shown on Fig. 
At the origin the three lowest eigenvalues are real and equal to 0, 32 and 60 which 
correspond to a free hyper-spheric real spectrum of K{K -\- 4) with K = 0,4, 6. 
The number K=2 is forbidden by the bosonic symmetry requirement. 

At large distances the lowest eigenvalue diverges quadratically according to 
due to a bound two-body state. The other eigenvalues return back to the 
free real spectrum at large distances with one of them taking the place of the 
diverged lowest eigenvalue. 

The energy spectrum for this model system is shown on Fig. |21 The points 
which form two lines represent the discretized continuum spectrum. The line 
which starts from E=0 corresponds to the states with genuine three-body asymp- 
totics where all three particles are at the continuum. Another line which starts 
at the two-body bound state energy E^'^^ = —6.76 MeV corresponds to the con- 
tinuum states where two of the particles are confined in the bound two-body 
state while the third one is in the continuum. 

Using four angular eigenvalues and pmax ~100 fm we find a bound three-body 
ground state at Eq = —37.22 MeV and a resonance at Ei = — 5.96 — 0.40i MeV in 
agreement with J^l- The independece of the calculated energies upon the scaling 
angle is illustrated in Tabled 

The radial function for the resonance is shown on Fig. |3J At shorter distances 
the function looks like a bound first excited state as indicated by the single node 
in the real part at shorter distances. At larger distances though, unlike the bound 
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Table 1. The ground state energy Eq and the resonance energy Ei for different scahng angles 
9 for the model system 






Eo, MeV 


El, MeV 


0.30 


-37.221 


-5.968 - 0.400i 


0.35 


-37.220 


-5.962 - 0.404i 


0.40 


-37.221 


-5.963 - 0.40H 




10 20 30 40 50 

p (fm) 

Figure 3. The dominating hyper-radial function /i of the lowest resonance in the model three- 
boson system. The scaling angle is 9=0.4. 

state functions, the resonance function oscillates with decreasing amphtude due 
to the complex scaling in accordance with 

3.2 Application to ^He 

In this section we apply the complex scaled adiabatic hyper-spheric method to 
the 2"*" resonance in the halo nucleus ^He within the three-body a-|-n-|-n model. 

We use the same interactions as in j21| . namely gaussian potentials fitted 
to the low energy scattering data and a phenomenological gaussian three-body 
force. This force is introduced to account for the additional polarisation of the 
alpha particle due to the presence of another neutron. It simulates the small 
additional attraction in a + n + n coming from excitations of the a-particle. The 
mass of the a-particle is assumed to be 4x939 MeV. 

The four lowest angular eigenvalues are shown on Fig. 0] for the scaling angle 
^=0.2. The eigenvalues start at p=0 from the real numbers 12 and 32 which 
correspond to the free real hyper-spheric spectrum of K{K + 4) with the hyper- 
angular momentum K=2 and 4. One of the eigenvalues then forms an attractive 
pocket where the resonance resides. This is largely due to n-a p- waves where 
the potential is attractive. The other two eigenvalues relate to n-a s-waves with 
strong repulsion which leads to the two repulsive eigenvalues contributing much 
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Figure 4. The four lowest complex eigenvalues for the a+n+n system, J'^=2"'". The scaling 
angle is 6—0.2. 
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Figure 5. The complex energy spectrum of ''He 2"*" states for the scaling angle 6=0.2. The res- 
onance at _B = 0.842 - 0.041i MeV is indicated by an arrow. The line Im{E) = - tan(2e')Re(_B) 
represents the exact rotated continuum spectrum. 



less to the resonance wave-function. At large distances the eigenvalues return 
back to the free real spectrum as there are no two-body bound states in this 
system and the scaling angle is too small for the two-body resonances to become 
bound. 

The complex scaled spectrum for 2+ states in ^He is shown on Fig. El No 
bound three-body states but one narrow resonance are present with these quan- 
tum numbers. With four adiabatic channels and pmax ~100 fm we obtain the 

resonance at E = 0.842 — 0.041i MeV which is reasonably close to the re- 
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Table 2. Resonance energy E for different scaling angles 6 for the a + n + n system 



e 


E, 


MeV 


0.15 


0.8419 


- 0.0405i 


0.20 


0.8417 


- 0.0404Z 


0.25 


0.8416 


- 0.0405Z 



0.15 




0.05 



-0.05 



I " " I " " I " " I " " I " " I " " I ' 
10 20 30 40 50 60 70 80 
p (fm) 



Figure 6. The dominating hyper-radial function /i of the Q-|-n-|-n system for the = 2^ 
resonance for the scaling angle 6=0.2. 



suit E = 0.845 — 0.046i MeV obtained in |21j with the complex energy method 
and pmax=180 fm. A similar result, E = 0.81 — 0.065i, has been obtained in 
jl6j . also with complex scaling method but with a different interaction model. 
The independence of the resonance energy on the scaling angle is illustrated in 
Table m 

The largest hyper-radial function of the resonance is shown on Fig. IHl In the 
inner region the hyper-radial function resembles a bound ground state function 
with no nodes and a considerable concentration in the pocket region. Again 
we see a fingerprint of a complex scaled resonance, namely the wave-function 
oscillating at large distances with ever decreasing amplitude according to (|29|) . 

3.3 Application to the three-a system 

The narrow 0"*" resonance in ^^C nucleus with E = 0.38 MeV and F = 8.3ibl.O eV 
was described quantitatively in JH] within a three-a model using the complex 
energy method. Here we repeat the calculation using the same parameters but 
instead with the complex scaling method. Unlike the first two numerical examples 
this one now includes a long range Coulomb potential which significantly alters 
the asymptotics of the wave- function |18j| compared to short range potentials. 
The a-a potential is taken as in ^Sl 

V{r) = (125MeVPi=o + 20MeVPi=2)e-'^'/(^-^^^'^)' 
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Complex scaling of the hyper-spheric coordinates... 
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Figure 7. The three lowest complex scaled angular eigenvalues for the three-a system. The 
scaling angle 6^=0.11. 



where Pi is the projection operator onto a state with relative orbital momentum I. 
The three-body force is Vs{p) = — 96.8MeV exp(— p^/ (3.9fm)^) where the nucleon 
mass m=939 MeV was used as a scale in the definition (|10llll|) of the hyper-radius 
P- 

The angular eigenvalues for this system are shown on Fig. [3 The lowest 
eigenvalue starts from zero corresponding to K = 0. It contains mostly repulsive 
s- waves and is always positive without any attractive pocket. The next eigenvalue 
starts from A=32 (K = 4) as the hyper-angular quantum number = 2 is 
forbidden by symmetry requirements. It contains mostly attractive p-waves and 
does form a relatively deep attractive pocket. 

At large p the real and imaginary parts of the eigenvalues grow linearly in- 
stead of converging to a real constant as is the case with short range potentials. 
This linear growth is due to the long range Coulomb potential jl8j which reap- 
pears in the effective potentials after the eigenvalue is divided by p'^. 

The energy spectrum is shown on Fig. [H] for the scaling angle 9=0.1 and 
Pmax=QO fm. There is a bound ground state located at £"0 = —6.79 MeV and 
two resonances, the very narrow one at Ei = 0.330 MeV and a broader resonance 
at E2 = 4.30 - 0.32i MeV. 

It is somewhat difficult to tell the first resonance from a discretized contin- 
uum state judging from the energy spectrum plot alone. However the hyper-radial 
wave- function shown on Fig. IHl leaves no doubts that this is indeed a resonance. 
Shown is the dominating function /2 which corresponds to the second eigenvalue 
on Fig. 13 where the attractive pocket is present. The function is again concen- 
trated at smaller distances and has one node in the real part corresponding to 
the first excited state. At larger distances it basically falls off without oscillations 
like as if it were a true bound state. This behavior results from the large effective 



30.18MeVe-'-'/(2-85fm)2 



(60) 
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T 

-4-2 2 

ReE (MeV) 

Figure 8. The complex energy spectrum of the three-a 0"*" states for the scaling angle 9=0.11. 
The resonance at £^ = 0.33 MeV is indicated by an arrow. 
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Figure 9. The dominating hyper-radial function /2(p) for the three-a system for the lowest 0"*" 
resonance and 0=0.11. 
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Complex scaling ol the hyper-spheric coordinates. 



Table 3. The bound state energy Eq and resonance energies Ei and E2 of the three lowest 
states in 3a system 






Eo, MeV 


El, MeV 


E2, MeV 


0.11 


-6.789 


0.3300 


4.295 - 0.323i 


0.12 


-6.793 


0.3283 


4.294 - 0.323i 


0.13 


-6.790 


0.3297 


4.296 - 0.325i 


0.14 


-6.790 


0.3296 


4.300 - 0.321i 
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Figure 10. The dominating hyper-radial function /i(p) for the three-a system for the second 
0^ resonance and 61=0.11. 



Coulomb barrier provided by the linearly growing angular eigenvalues as shown 
on Fig. [71 

Another criterion to distinguish the resonances from the continuum states is 
that when the scaling angle 6 is changed the whole line of continuum states ro- 
tates according to (|31|) while the resonances (and bound states) are independent 
of the scaling angle as illustrated in Table 01 

The energies obtained in this calculation compare well to the values of - 
6.81 MeV for the ground state and 380-0. 02i keV for the first resonance obtained 
in JHI with the complex-energy method. The source of the difference is that in 
[TH] only two of the lowest angular eigenvalues were used while here we use 
as many eigenvalues as needed for convergence - typically 4-5. Again for these 
numerical illustrations we have more or less arbitrarily chosen about 100 mesh 
points in the hyperradius for all examples. No fine tuning has been made in order 
to achieve higher accuracy. The present numerical error is estimated to be about 
10 keV and, therefore, the width of the order of 0.01 keV of the first extremely 
narrow resonance is beyond the capabilities of an untuned mesh. 

Another resonance, which was not calculated in 18 , is clearly seen in the 
spectrum with the energy E = 4.30 — 0.32i MeV, that is with the excitation 
energy £'*=11 MeV and the width 1^=0.6 MeV. It apparently corresponds to the 
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experimentally known broader second excited 0^ state with £'*=10.3ib0.3 MeV 
and r=3±0.7 MeV (22^. 

The dominating radial function, which for this resonance corresponds to the 
first eigenvalue on Fig. [7| is shown on Fig. IIUI The real part of the wavefunction 
has no nodes in the pocket region. Therefore this resonance represents the lowest 
state in the hyper-angular eigenvalue number 1. This eigenvalue also forms a 
pocket, see Fig. although not as deep as that of eigenvalue number 2. 

4 Conclusion 

We have implemented complex scaling of coordinates within the hyper-spheric 
adiabatic approach to Faddeev equations. We have derived the necessary equa- 
tions and investigated the asymptotic spectrum of the hyper-angular eigenvalues. 

The complex-scaled hyper-radial wave-functions fall off exponentially at large 
distances which allows the simple boundary condition of vanishing functions to be 
used in the search for the resonances. Alternatively, it is possible to use square 
integrable functions as a basis set, thus avoiding the less effective shoot-and- 
match method. Using such a basis expansion has an additional advantage that 
all bound states and resonances are calculated in one diagonalization run while 
otherwise one has to search separately for every single pole of the S-matrix. Again 
with basis expansion in addition to resonances and bound states one also gets in 
the same run the (discretized) continuum wave functions. Otherwise continuum 
wave functions have to be calculated separately. 

The simplified boundary conditions are especially important for the three 
positively charged particles where otherwise the non-scaled boundary condition 
involves complicated functions with slow asymptotic convergence. 

The drawbacks of the complex scaling method is first that it requires the 
numerically slower complex arithmetics to be used and second that the ranges of 
short-range potentials effectively become larger after complex scaling with the 
numerically stronger demand of more accurate treatment at larger distances. 

We have illustrated the viability of the method by applications to three sys- 
tems, one benchmark model system and two realistic systems, one with short- 
range potentials and the other in addition also including the long-range Coulomb 
potentials. In conclusion the complex-scaled hyper-spheric approach to Faddeev 
equations is a viable alternative method for calculating three-body resonances 
and continuum spectrum states. 
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